Spectral, Compositional, and Physical Properties of the Upper Murray Formation and Vera Rubin Ridge, Gale Crater, Mars.

During 2018 and 2019, the Mars Science Laboratory Curiosity rover investigated the chemistry, morphology, and stratigraphy of Vera Rubin ridge (VRR). Using orbital data from the Compact Reconnaissance Imaging Spectrometer for Mars, scientists attributed the strong 860 nm signal associated with VRR to the presence of red crystalline hematite. However, Mastcam multispectral data and CheMin X-ray diffraction (XRD) measurements show that the depth of the 860 nm absorption is negatively correlated with the abundance of red crystalline hematite, suggesting that other mineralogical or physical parameters are also controlling the 860 nm absorption. Here, we examine Mastcam and ChemCam passive reflectance spectra from VRR and other locations to link the depth, position, and presence or absence of iron-related mineralogic absorption features to the XRD-derived rock mineralogy. Correlating CheMin mineralogy to spectral parameters showed that the ~860 nm absorption has a strong positive correlation with the abundance of ferric phyllosilicates. New laboratory reflectance measurements of powdered mineral mixtures can reproduce trends found in Gale crater. We hypothesize that variations in the 860 nm absorption feature in Mastcam and ChemCam observations of VRR materials are a result of three factors: (1) variations in ferric phyllosilicate abundance due to its ~800-1,000 nm absorption; (2) variations in clinopyroxene abundance because of its band maximum at ~860 nm; and (3) the presence of red crystalline hematite because of its absorption centered at 860 nm. We also show that relatively small changes in Ca-sulfate abundance is one potential cause of the erosional resistance and geomorphic expression of VRR.

Structure and dynamics of water and lipid molecules in charged anionic DMPG lipid bilayer membranes.

Molecular dynamics simulations have been used to investigate the influence of the valency of counter-ions on the structure of freestanding bilayer membranes of the anionic 1,2-dimyristoyl-sn-glycero-3-phosphoglycerol (DMPG) lipid at 310 K and 1 atm. At this temperature, the membrane is in the fluid phase with a monovalent counter-ion and in the gel phase with a divalent counter-ion. The diffusion constant of water as a function of its depth in the membrane has been determined from mean-square-displacement calculations. Also, calculated incoherent quasielastic neutron scattering functions have been compared to experimental results and used to determine an average diffusion constant for all water molecules in the system. On extrapolating the diffusion constants inferred experimentally to a temperature of 310 K, reasonable agreement with the simulations is obtained. However, the experiments do not have the sensitivity to confirm the diffusion of a small component of water bound to the lipids as found in the simulations. In addition, the orientation of the dipole moment of the water molecules has been determined as a function of their depth in the membrane. Previous indirect estimates of the electrostatic potential within phospholipid membranes imply an enormous electric field of 10(8)-10(9) V m(-1), which is likely to have great significance in controlling the conformation of translocating membrane proteins and in the transfer of ions and molecules across the membrane. We have calculated the membrane potential for DMPG bilayers and found ∼1 V (∼2 ⋅ 10(8) V m(-1)) when in the fluid phase with a monovalent counter-ion and ∼1.4 V (∼2.8 ⋅ 10(8) V m(-1)) when in the gel phase with a divalent counter-ion. The number of water molecules for a fully hydrated DMPG membrane has been estimated to be 9.7 molecules per lipid in the gel phase and 17.5 molecules in the fluid phase, considerably smaller than inferred experimentally for 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) membranes but comparable to the number inferred for 1,2-dilauroyl-sn-glycero-3-phosphoethanolamine (DLPE) membranes. Some of the properties of the DMPG membrane are compared with those of the neutral zwitterionic DMPC bilayer membrane at 303 K and 1 atm, which is the same reduced temperature with respect to the gel-to-fluid transition temperature as 310 K is for the DMPG bilayer membrane.

Diffusion of water and selected atoms in DMPC lipid bilayer membranes.

Molecular dynamics simulations have been used to determine the diffusion of water molecules as a function of their position in a fully hydrated freestanding 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) bilayer membrane at 303 K and 1 atm. The diffusion rate of water in a ∼10 Šthick layer just outside the membrane surface is reduced on average by a factor of ∼2 relative to bulk. For water molecules penetrating deeper into the membrane, there is an increasing reduction in the average diffusion rate with up to one order of magnitude decrease for those deepest in the membrane. A comparison with the diffusion rate of selected atoms in the lipid molecules shows that ∼6 water molecules per lipid molecule move on the same time scale as the lipids and may therefore be considered to be tightly bound to them. The quasielastic neutron scattering functions for water and selected atoms in the lipid molecule have been simulated and compared to observed quasielastic neutron scattering spectra from single-supported bilayer DMPC membranes.

Steric hindrance as a basis for structure-based design of selective inhibitors of protein-tyrosine phosphatases.

Utilizing structure-based design, we have previously demonstrated that it is possible to obtain selective inhibitors of protein-tyrosine phosphatase 1B (PTP1B). A basic nitrogen was introduced into a general PTP inhibitor to form a salt bridge to Asp48 in PTP1B and simultaneously cause repulsion in PTPs containing an asparagine in the equivalent position [Iversen, L. F., et al. (2000) J. Biol. Chem. 275, 10300-10307]. Further, we have recently demonstrated that Gly259 in PTP1B forms the bottom of a gateway that allows easy access to the active site for a broad range of substrates, while bulky residues in the same position in other PTPs cause steric hindrance and reduced substrate recognition capacity [Peters, G. H., et al. (2000) J. Biol. Chem. 275, 18201-18209]. The current study was undertaken to investigate the feasibility of structure-based design, utilizing these differences in accessibility to the active site among various PTPs. We show that a general, low-molecular weight PTP inhibitor can be developed into a highly selective inhibitor for PTP1B and TC-PTP by introducing a substituent, which is designed to address the region around residues 258 and 259. Detailed enzyme kinetic analysis with a set of wild-type and mutant PTPs, X-ray protein crystallography, and molecular modeling studies confirmed that selectivity for PTP1B and TC-PTP was achieved due to steric hindrance imposed by bulky position 259 residues in other PTPs.

Influence of a lipid interface on protein dynamics in a fungal lipase.

Lipases catalyze lipolytic reactions and for optimal activity they require a lipid interface. To study the effect of a lipid aggregate on the behavior of the enzyme at the interfacial plane and how the aggregate influences an attached substrate or product molecule in time and space, we have performed molecular dynamics simulations. The simulations were performed over 1 to 2 ns using explicit SPC water. The interaction energies between protein and lipid are mainly due to van der Waals contributions reflecting the hydrophobic nature of the lipid molecules. Estimations of the protonation state of titratable residues indicated that the negative charge on the fatty acid is stabilized by interactions with the titratable residues Tyr-28, His-143, and His-257. In the presence of a lipid patch, the active site lid opens wider than observed in the corresponding simulations in an aqueous environment. In that lid conformation, the hydrophobic residues Ile-85, Ile-89, and Leu-92 are embedded in the lipid patch. The behavior of the substrate or product molecule is sensitive to the environment. Entering and leaving of substrate molecules could be observed in presence of the lipid patch, whereas the product forms strong hydrogen bonds with Ser-82, Ser-144, and Trp-88, suggesting that the formation of hydrogen bonds may be an important contribution to the mechanism by which product inhibition might take place.

Dynamics of the substrate binding pocket in the presence of an inhibitor covalently attached to a fungal lipase.

To gain insight into the mobility of the occupied ligand-binding pocket of the Rhizomucor miehei lipase we have conducted a rigorous molecular dynamics analysis. The covalently attached inhibitor, ethylhexylphosphonate, was employed as a mimic of the putative tetrahedral intermediate in the esterolytic reaction. Our results show that in this lipase, ligand recognition is influenced by the flexibility of the binding pocket, a feature that is common to many other enzymes. Several regions around the active site were found to move significantly to adapt to the inhibitor. These motions are correlated to the flexibility of the inhibitor. In particular, the hexyl chain of the inhibitor shows considerable mobility, and adjacent residues in the binding cleft accommodate to this flexibility. Pronounced fluctuations in the binding pocket induced by the flexibility of the inhibitor are observed in the hinge region F79-S82, the active site loop region W88-V95 and the protein regions P209-F215/H257-Y260. The flexibility in the regions F79-S82 and H257-Y260, where the shorter ethyl chain is located, indicates that additional space in this binding cleft region is available for accommodating a larger moiety. Fluctuations in the region W88-V95 and P209-F215 are due to the relatively short flexible hexyl carbon chain. This part of the binding pocket could be stiffened by the presence of a longer carbon chain. Though the inhibitor is covalently attached through the phosphonate moiety, interaction of the remainder of the molecule and the enzyme are determined by hydrophobic interactions, where the Van der Waals energies are approximately 25% lower than the electrostatic contributions.

Dipolar and chain-linking effects on the rheology of grafted chains in a nanopore under shear at different grafting densities.

Nonequilibrium molecular dynamics simulations are applied to investigate the rheological properties of coplanar nanopore systems of amphiphilic chain molecules with the tails grafted to the walls of the nanopore and with the head-group ends immersed in a solvent inside the nanopore. In particular, the effects of modifying the interaction between the amphiphilic head-groups by repulsive dipolar interactions or directly covalently linking pairs of chains at the head-groups have been studied. Different grafting densities are considered. The chains are modeled by a harmonic bead-spring model, and all particles interact through the repulsive part of a shifted Lennard-Jones potential. Head-group linking is also governed by a bead-spring potential. A harmonic potential models the lattice vibrations of the atomic boundaries. The rheological properties are studied by a shearing process in which the heat generated is conducted away from the system through the walls by applying a Nosé-Hoover thermostat. Computed geometric parameters such as average chain length and average tilt angle indicate reduction in chain flexibility at large dipole moments. Dipolar repulsion is found to broaden the density profiles of the solvent. This effect is opposed by chain linking. For increasing head-group repulsion, the amphiphile-solvent interfaces become less diffusive that leads to systematic variations in viscosities with increasing dipole moments. Friction forces become stronger at large grafting density and for larger dipole moments. The changes in rheological properties for fixed grafting density are essentially governed by the change in the response of the normal pressure to the applied shear field. The velocity gradients depend strongly on the degree of grafting density but appear to be less sensitive to the strength of the interactions between the head groups.

Structure-based design of a low molecular weight, nonphosphorus, nonpeptide, and highly selective inhibitor of protein-tyrosine phosphatase 1B.

Several protein-tyrosine phosphatases (PTPs) have been proposed to act as negative regulators of insulin signaling. Recent studies have shown increased insulin sensitivity and resistance to obesity in PTP1B knockout mice, thus pointing to this enzyme as a potential drug target in diabetes. Structure-based design, guided by PTP mutants and x-ray protein crystallography, was used to optimize a relatively weak, nonphosphorus, nonpeptide general PTP inhibitor (2-(oxalyl-amino)-benzoic acid) into a highly selective PTP1B inhibitor. This was achieved by addressing residue 48 as a selectivity determining residue. By introducing a basic nitrogen in the core structure of the inhibitor, a salt bridge was formed to Asp-48 in PTP1B. In contrast, the basic nitrogen causes repulsion in other PTPs containing an asparagine in the equivalent position resulting in a remarkable selectivity for PTP1B. Importantly, this was accomplished while retaining the molecular weight of the inhibitor below 300 g/mol.

Residue 259 is a key determinant of substrate specificity of protein-tyrosine phosphatases 1B and alpha.

The aim of this study was to define the structural elements that determine the differences in substrate recognition capacity of two protein-tyrosine phosphatases (PTPs), PTP1B and PTPalpha, both suggested to be negative regulators of insulin signaling. Since the Ac-DADE(pY)L-NH(2) peptide is well recognized by PTP1B, but less efficiently by PTPalpha, it was chosen as a tool for these analyses. Calpha regiovariation analyses and primary sequence alignments indicate that residues 47, 48, 258, and 259 (PTP1B numbering) define a selectivity-determining region. By analyzing a set of DADE(pY)L analogs with a series of PTP mutants in which these four residues were exchanged between PTP1B and PTPalpha, either in combination or alone, we here demonstrate that the key selectivity-determining residue is 259. In PTPalpha, this residue is a glutamine causing steric hindrance and in PTP1B a glycine allowing broad substrate recognition. Significantly, replacing Gln(259) with a glycine almost turns PTPalpha into a PTP1B-like enzyme. By using a novel set of PTP inhibitors and x-ray crystallography, we further provide evidence that Gln(259) in PTPalpha plays a dual role leading to restricted substrate recognition (directly via steric hindrance) and reduced catalytic activity (indirectly via Gln(262)). Both effects may indicate that PTPalpha regulates highly selective signal transduction processes.

Molecular dynamics simulations of protein-tyrosine phosphatase 1B. II. substrate-enzyme interactions and dynamics.

Molecular dynamics simulations of protein tyrosine phosphatase 1B (PTP1B) complexed with the phosphorylated peptide substrate DADEpYL and the free substrate have been conducted to investigate 1) the physical forces involved in substrate-protein interactions, 2) the importance of enzyme and substrate flexibility for binding, 3) the electrostatic properties of the enzyme, and 4) the contribution from solvation. The simulations were performed for 1 ns, using explicit water molecules. The last 700 ps of the trajectories was used for analysis determining enthalpic and entropic contributions to substrate binding. Based on essential dynamics analysis of the PTP1B/DADEpYL trajectory, it is shown that internal motions in the binding pocket occur in a subspace of only a few degrees of freedom. In particular, relatively large flexibilities are observed along several eigenvectors in the segments: Arg(24)-Ser(28), Pro(38)-Arg(47), and Glu(115)-Gly(117). These motions are correlated to the C- and N-terminal motions of the substrate. Relatively small fluctuations are observed in the region of the consensus active site motif (H/V)CX(5)R(S/T) and in the region of the WPD loop, which contains the general acid for catalysis. Analysis of the individual enzyme-substrate interaction energies revealed that mainly electrostatic forces contribute to binding. Indeed, calculation of the electrostatic field of the enzyme reveals that only the field surrounding the binding pocket is positive, while the remaining protein surface is characterized by a predominantly negative electrostatic field. This positive electrostatic field attracts negatively charged substrates and could explain the experimentally observed preference of PTP1B for negatively charged substrates like the DADEpYL peptide.

Computational analysis of chain flexibility and fluctuations in Rhizomucor miehei lipase.

We have performed molecular dynamics simulation of Rhizomucor miehei lipase (Rml) with explicit water molecules present. The simulation was carried out in periodic boundary conditions and conducted for 1. 2 ns in order to determine the concerted protein dynamics and to examine how well the essential motions are preserved along the trajectory. Protein motions are extracted by means of the essential dynamics analysis method for different lengths of the trajectory. Motions described by eigenvector 1 converge after approximately 200 ps and only small changes are observed with increasing simulation time. Protein dynamics along eigenvectors with larger indices, however, change with simulation time and generally, with increasing eigenvector index, longer simulation times are required for observing similar protein motions (along a particular eigenvector). Several regions in the protein show relatively large fluctuations and in particular motions in the active site lid and the segments Thr57-Asn63 and the active site hinge region Pro101-Gly104 are seen along several eigenvectors. These motions are generally associated with glycine residues, while no direct correlations are observed between these fluctuations and the positioning of prolines in the protein structure. The partial opening/closing of the lid is an example of induced fit mechanisms seen in other enzymes and could be a general mechanism for the activation of Rml.

Molecular dynamics simulations of protein-tyrosine phosphatase 1B. I. ligand-induced changes in the protein motions.

Activity of enzymes, such as protein tyrosine phosphatases (PTPs), is often associated with structural changes in the enzyme, resulting in selective and stereospecific reactions with the substrate. To investigate the effect of a substrate on the motions occurring in PTPs, we have performed molecular dynamics simulations of PTP1B and PTP1B complexed with a high-affinity peptide DADEpYL, where pY stands for phosphorylated tyrosine. The peptide sequence is derived from the epidermal growth factor receptor (EGFR988-993). Simulations were performed in water for 1 ns, and the concerted motions in the protein were analyzed using the essential dynamics technique. Our results indicate that the predominately internal motions in PTP1B occur in a subspace of only a few degrees of freedom. Upon substrate binding, the flexibility of the protein is reduced by approximately 10%. The largest effect is found in the protein region, where the N-terminal of the substrate is located, and in the loop region Val198-Gly209. Displacements in the latter loop are associated with the motions in the WPD loop, which contains a catalytically important aspartic acid. Estimation of the pKa of the active-site cysteine along the trajectory indicates that structural inhomogeneity causes the pKa to vary by approximately +/-1 pKa unit. In agreement with experimental observations, the active-site cysteine is negatively charged at physiological pH.

Analysis of the dynamics of rhizomucor miehei lipase at different temperatures.

The dynamics of Rhizomucor miehei lipase has been studied by molecular dynamics simulations at temperatures ranging from 200-500K. Simulations carried out in periodic boundary conditions and using explicit water molecules were performed for 400 ps at each temperature. Our results indicate that conformational changes and internal motions in the protein are significantly influenced by the temperature increase. With increasing temperature, the number of internal hydrogen bonds decreases, while surface accessibility, radius of gyration and the number of residues in random coil conformation increase. In the temperature range studied, the motions can be described in a low dimensional subspace, whose dimensionality decreases with increasing temperature. Approximately 80% of the total motion is described by the first (i) 80 eigenvectors at T=200K, (ii) 30 eigenvectors at T=300K and (iii) 10 eigenvectors at T=400K. At high temperature, the alpha-helix covering the active site in the native Rhizomucor miehei lipase, the helix at which end the active site is located, and in particular, the loop (Gly35-Lys50) show extensive flexibility.

Assignment of side-chain conformation using adiabatic energy mapping, free energy perturbation, and molecular dynamic simulations.

NMR spectroscopic analysis of the C-terminal Kunitz domain fragment (alpha3(VI)) from the human alpha3-chain of type VI collagen has revealed that the side chain of Trp21 exists in two unequally populated conformations. The major conformation (M) is identical to the conformation observed in the X-ray crystallographic structure, while the minor conformation (m) cannot structurally be resolved in detail by NMR due to insufficient NOE data. In the present study, we have applied: (1) rigid and adiabatic mapping, (2) free energy simulations, and (3) molecular dynamic simulations to elucidate the structure of the m conformer and to provide a possible pathway of the Trp21 side chain between the two conformers. Adiabatic energy mapping of conformations of the Trp21 side chain obtained by energy minimization identified two energy minima: One corresponding to the conformation of Trp21 observed in the X-ray crystallographic structure and solution structure of alpha3(VI) (the M conformation) and the second corresponding to the m conformation predicted by NMR spectroscopy. A transition pathway between the M and m conformation is suggested. The free-energy difference between the two conformers obtained by the thermodynamic integration method is calculated to 1.77+/-0.7 kcal/mol in favor of the M form, which is in good agreement with NMR results. Structural and dynamic properties of the major and minor conformers of the alpha3(VI) molecule were investigated by molecular dynamic. Essential dynamics analysis of the two resulting 800 ps trajectories reveals that when going from the M to the m conformation only small, localized changes in the protein structure are induced. However, notable differences are observed in the mobility of the binding loop (residues Thr13-Ile18), which is more flexible in the m conformation than in the M conformation. This suggests that the reorientation of Trp2 might influence the inhibitory activity against trypsin, despite the relative large distance between the binding loop and Trp21.

Active serine involved in the stabilization of the active site loop in the Humicola lanuginosa lipase.

We have investigated the binding properties of and dynamics in Humicola lanuginosa lipase (Hll) and the inactive mutant S146A (active Ser146 substituted with Ala) using fluorescence spectroscopy and molecular dynamics simulations, respectively. Hll and S146A show significantly different binding behavior for phosphatidylcholine (PC) and phosphatidylglycerol (PG) liposomes. Generally, higher binding affinity is observed for Hll than the S146A mutant. Furthermore, depending on the matrix, the addition of the transition state analogue benzene boronic acid increases the binding affinity of S146A, whereas only small changes are observed for Hll suggesting that the active site lid in the latter opens more easily and hence more lipase molecules are bound to the liposomes. These observations are in agreement with molecular dynamics simulations and subsequent essential dynamics analyses. The results reveal that the hinges of the active site lid are more flexible in the wild-type Hll than in S146A. In contrast, larger fluctuations are observed in the middle region of the active site loop in S146A than in Hll. These findings reveal that the single mutation (S146A) of the active site serine leads to substantial conformational alterations in the H. lanuginosa lipase and different binding affinities.

Electrostatic evaluation of the signature motif (H/V)CX5R(S/T) in protein-tyrosine phosphatases.

The catalytic activity of protein-tyrosine phosphatases (PTPs) is mediated by a cysteine side chain which carries out a nucleophilic attack initiating the phosphate cleavage. Experimentally, it has been observed that the active site cysteine has a remarkably low pKa. In the present study, we have investigated the origin of the low pKa by analyzing the electrostatic properties of four different protein-tyrosine phosphatases: Yersinia PTP (bacteria), PTP1B (human), VHR (human), and low molecular weight phosphatase (bovine). These phosphatases have very low sequence homology and show very low structural similarity. However, they share a common active site motif [the (H/V)CX5R(S/T) sequence] which adopts a unique loop structure. We have applied the so-called single site titration method, which is based on the Poisson-Boltzmann methodology, to (i) study the influence of the architecture of the (H/V)CX5R(S/T) loop on the pKa of the active cysteine and (ii) examine which parts of the active site region stabilize the ionized form of the cysteine. Our results indicate that the architecture of the (H/V)CX5R(S/T) loop has a major impact on the low pKa of the active cysteines. The orientation of the microdipoles generated by the partial charges of the backbone atoms (i.e., the CONHCalpha atoms) is essential for maintaining the low pKa. Further, the electrostatic field generated by these microdipoles has a larger impact than the electrostatic dipole generated by the central alpha-helix. Interactions of the active cysteine with other ionizable side chains play a minor role in stabilizing the thiolate anion. The only ionizable side chain significantly influencing the pKa of the active site cysteine is the arginine, which is an important part of the consensus sequence.

Computational studies of the activation of lipases and the effect of a hydrophobic environment.

We have investigated the activation pathway of three wild type lipases and three mutants using molecular dynamics techniques combined with a constrained mechanical protocol. The activation of these lipases involves a rigid body hinge-type motion of a single helix, which is displaced during activation to expose the active site and give access to the substrate. Our results suggest that the activation of lipases is enhanced in a hydrophobic environment as is generally observed in experiments. The energy gain upon activation varies between the different lipases and depends strongly on the distribution of the charged residues in the activating loop region. In a low dielectric constant medium (such as a lipid environment), the electrostatic interactions between the residues located in the vicinity of the activating loop (lipid contact zone) are dominant and determine the activation of the lipases. Calculations of the pKas qualitatively indicate that some titratable residues experience significant pK shifts upon activation. These calculations may provide sufficient details for an understanding of the origin and magnitude of a given electrostatic effect and may provide an avenue for exploring the activation pathway of lipases.

Essential dynamics of lipase binding sites: the effect of inhibitors of different chain length.

The biochemical activity of enzymes, such as lipases, is often associated with structural changes in the enzyme resulting in selective and stereospecific reactions with the substrate. To investigate the effect of a substrate and its chain length on the dynamics of the enzyme, we have performed molecular dynamics simulations of the native Rhizomucor miehei lipase (Rml) and lipase-dialkylphosophate complexes, where the length of the alkyl chain ranges from two to 10 carbon atoms. Simulations were performed in water and trajectories of 400 ps were used to analyse the essential motions in these systems. Our results indicate that the internal motions of the Rml and Rml complexes occur in a subspace of only a few degrees of freedom. A high flexibility is observed in solvent-exposed segments, which connect beta-sheets and helices. In particular, loop regions Gly35-Lys50 and Thr57-Asn63 fluctuate extensively in the native enzyme. Upon activation and binding of the inhibitor, involving the displacement of the active site loop, these motions are considerably suppressed. With increasing chain length of the inhibitor, the fluctuations in the essential subspace increase, levelling off at a chain length of 10, which corresponds to the size of the active-site groove.

Dynamics of proteins in different solvent systems: analysis of essential motion in lipases.

We have investigated the effect of different solvents on the dynamics of Rhizomucor miehei lipase. Molecular dynamics simulations were performed in water, methyl hexanoate, and cyclohexane. Analysis of the 400-ps trajectories showed that the solvent has a pronounced effect on the geometrical properties of the protein. The radius of gyration and total accessibility surface decrease in organic solvents, whereas the number of hydrogen bonds increases. The essential motions of the protein in different solvents can be described in a low-dimensional "essential subspace," and the dynamic behavior in this subspace correlates with the polarity of the solvent. Methyl hexanoate, which is a substrate for R. miehei lipase, significantly increases the fluctuations in the active-site loop. During the simulation, a methyl hexanoate entered the active-site groove. This observation provides insight into the possible docking mechanism of the substrate.